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The diverse phenomena of cavity nonlinear 
optics [lj provide a rich basis for fundamental 
studies of dissipative nonlinear dynamics [2l [3] 
and for the design of photonic signal process- 
ing devices |4j. Recent experiments exploit- 
ing resonant atomic nonlinearities [5HT0J and/or 
nanophotonic cavities |HH17j have demonstrated 
that such phenomena can occur at very low 
(femtojoule attojoule) energy scales in prototyp- 
ical systems, reaching down to the quantum- 
physical few photon regime [18, 19j and raising 
intriguing prospects for corresponding ultra-low 
power photonic information technology |20H22| . 
One of the less-explored dimensions of this de- 
veloping scenario is the surprising complexity of 
dynamical behaviors, beyond simple thresholding 
and bistability, that can be achieved at low energy 
scales in cavities incorporating two-level atoms or 
comparable solid-state emitters. In this Letter 
we present experimental data on sub-critical Hopf 
bifurcation phenomena in an atom-cavity system 
and demonstrate a novel device configuration that 
latches a logic level and provides rf modulation of 
the output signal at the same time. 

The general paradigm of two-level emitters coupled to 
an isolated optical mode can be realized in a wide range 
of gas, liquid, and solid-state systems with various types 
of optical cavities; our experiment uses a 12mm Fabry- 
Perot resonator in ultrahigh vacuum with an ensemble 
of laser-cooled 133 Cs atoms placed between the mirrors. 
The atoms are collected and cooled on a chip, and then 
transported to the cavity magnetically [23] (see Meth- 
ods section for detailed experimental procedures). We 
utilize only the |F=niF=4) to |F'=m' F =5) transition of 
the D2 line at 852nm where the atoms are well modeled 
as two- level systems. The maximum atom-cavity cou- 
pling constant, the atomic polarization decay rate, and 
the measured cavity field decay rate are (g , 7_i_, k)/2tt = 
(0.72, 2.6, 2.3) MHz, respectively. 

We model the system composed of a cavity and N 
weakly-coupled, noninteracting atoms by the following 



semiclassical equations of motion [53] : 

N 

(a) = -( K + iA c )(a)+Y / 9j(<7 i -)+£(t) (la) 

3=1 

<*L) = -(7± + *A A ) (ol) + 9j (a) (ai) (lb) 
(oi) = -2 7 x«4) + 1) - %Re{ (<r J + ) (a) } (lc) 

Here, a is the intracavity field annihilation operator, ct± 
are the rasing/lowering operators for the j-th atom, and 
o-j, = [er J + ,er J _]. The rates gj,jj_,n are the atom-cavity 
coupling strength for the j-th atom, the atomic polar- 
ization decay rate, and the cavity field decay rate, re- 
spectively. Ac and Aa are the cavity and the atom de- 
tunings relative to the drive (laser) frequency, defined 
as A c = w C av - ^ las and A A = w atom - cj las , and £(t) 
is the drive field amplitude. The above Maxwell-Bloch- 
type equations are derived from an unconditional Jaynes- 
Cummings master equation in the same manner as in [19] 
(see also [2F]), but with a slight generalization that in- 
corporates distinct coupling strengths for the individual 
atoms. We note that the 'factorization' approximations 
made in the derivation, such as (aer^_) w (a) (o^_), mean 
that entanglement between the atoms and the field is ig- 
nored. This is justified because the coupling strengths for 
individual atoms are small, while the combined coupling 
strength as an ensemble is still large enough to cause a 
strong nonlinearity. 

In our experiment we assume a large number of atoms 
distributed uniformly inside the cavity, so that the val- 
ues of gj follow directly from the gaussian standing-wave 
profile of the resonant field mode (TEMoo). Under this 
assumption the atomic ensemble is characterized by a 
single parameter N c g = J2f=i(9j/9o) 2 , where go is the 
maximum coupling constant at the anti-node center of 
the cavity. This simplification applies best for a cavity 
filled with a homogeneous dispersive/absorptive medium, 
a scenario for which extensive studies have already been 
performed using the Maxwell-Bloch Equations (MBEs) 
[2 [2S] . The method used here differs from the conven- 
tional treatment of the MBEs in that we do not simplify 



the equations (lal-(lc) any further by the use of 'col- 
lective' atomic variables. We instead keep track of the 
evolution of ~ 10 3 individual atomic variables for better 
accuracy, at the cost of increased computational load. 



By setting the time derivatives in equations (la|-(lc 
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to zero one can solve for the steady state values of (a), 
(<r J _), and ((T 3 Z ). The independent variable £ is then ex- 
pressed as a function of (a) to derive the input-output 
characteristic of the system, and the stability of any (hy- 
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FIG. 1. Steady state nonlinear input-output character- 
istics, a, Threshold behavior of the system with N a g = 130, 
Ac = Aa = 0. The black line is the stable steady state solu- 
tion derived from the model, and red (green) line is a typical 
data trace from the measurement when the input power is 
gradually increased (decreased), b, Hysteretic behavior with 
TVpff = 470, Ac = Aa = 2n ■ 8MHz. Color scheme is the same 
as a except for the new black dotted line that indicates the un- 
stable steady state solution, c, Experimental demonstration 
of latching using bistable system parameters (as in b). With 
the input power constant at 14.5nW the system maintains its 
current state. When a pulse is applied the system state is set 
high or low according to the pulse type (positive/negative). 



perbolic) equilibrium point can be determined by local 
linearization. Alternatively, equations can be numer- 
ically integrated to produce sample trajectories that re- 
flect nonlinear dynamical phenomena such as limit cycles, 
chaos ([23 EH]), etc. In this Letter we make use of both 
equilibrium and transient data to assess the accuracy of 
our theoretical modeling approach. 

The accuracy of our model in a regime of relatively sim- 
ple dynamic behavior is illustrated in Figs. [l^,b, which 
display the measured equilibrium input-output response 
of our atom-cavity system for two different parameter 
sets {N c g, Ac,Aa} together with corresponding steady 
state solutions of equations For the parameters of 
Fig. [TJi (iV cff = 130, A c = A A = 0) a soft thresh- 
olding behavior is observed with a corner near 6 nW 
input power. With increased atom number and a de- 
tuned probe (N cS = 470, A c = A A = 2?r • 8MHz) as in 
Fig. we observe classic bistability and hysteresis with 
coexisting stable equilibria from roughly 14 to 15.5 nW 
input power. For these measurements the input power 
was swept slowly (20kHz) compared to the dynamical 
rates of the system (g; , 7j_ , k) to probe the adiabatic re- 
sponse. Figures fl&,b each show an individual (but pre- 
cisely repeatable) data trace of ~ 50/is duration that was 
recorded continuously as the input power was swept up 
and down. 



The theory curves in Figs. [T^,b were computed us- 
ing manually-adjusted values of 7V e ff and £(t) to improve 
agreement with our data, but the fitted values differ from 
the independently measured ones by no more than 10% 
(for the measurement of N c g, refer to [!§]). This pro- 
cedure and statement of calibration accuracy apply to 
all remaining data-theory comparison plots in this Let- 
ter as well. Some discrepancy between the model and 
the measurement is expected, mainly because the mag- 
netic trap from which the atoms are released just before 
the measurement has a similar size to the cavity mode 
(waist wq = 66/im), and therefore it is difficult to ensure 
that the atoms are truly uniformly distributed when the 
measurement is taken. Up to these ~ 10% inaccuracies 
however our model has proven to be a reliable predictor 
of the parameter values at which interesting nonlinear 
dynamical phenomena should be observed. 

After confirming the parameters for simple bistable 
input-output behavior in our atom-cavity system, we per- 
formed an initial proof-of-principle experiment to demon- 
strate an all-optical set-reset (SR) latch, which could for 
example be used as an optical memory bit. Typical data 
is shown in Fig. [TJ:. Here we envision a device configura- 
tion in which the set /reset control beam shares an opti- 
cal input channel with the bias power beam (for example 
they could be combined by a beam splitter) so that they 
interfere constructively or destructively depending on the 
phase of the control beam. With this picture in mind, we 
demonstrated the set /reset control by superposing posi- 
tive and negative pulses onto the laser beam that drives 
the atom-cavity system (upper blue trace of Fig. [TJ;) . A 
positive pulse switches the response to the upper branch 
of the hysteresis curve, where it latches until a negative 
pulse resets it back to the lower branch. Our system re- 
quired 15nW optical bias power and the energy used for 
each control pulse was roughly 5nW x 0.2/is = If J. In 
terms of the switching energy, this system is better than 
the nano-cavity based all-optical bistable devices devel- 
oped so far [T3], but we also find that this energy scale 
is still much larger than that of the optical shot noise 
(quantum fluctuations), which can destabilize the device 
and make it unreliable [TU] . 

With an appropriate choice of external parameters 
Nes, Ac and Aa guided by theory, we also probed the 
atom-cavity dynamics in an unstable regime and ob- 
served "self-pulsing" behavior 2, 30-32 in which the 
global attractor is a limit cycle. In Fig. [2^, we show the 
system making transitions between a stable equilibrium 
point and a limit cycle as we sweep the driving input 
power. When the input power is large, the power of the 
output beam oscillates significantly at a frequency in the 
range of - 21 - 24 MHz. 

Our model predicts a range of input powers (~ 1150 — 
1200nW) for which a stable equilibrium point and a sta- 
ble limit cycle coexist. Starting from low input, as the 
system is driven past the sub-critical Hopf bifurcation 
point pi (~1200nW, shown in Figs.[2ji,e) the equilibrium 
point becomes unstable and the system jumps to the the 
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FIG. 2. Hopf bifurcations and limit cycles. Experimen- 
tal parameters 7V cft = 2100, Ac/2tt = -20MHz, A a /2tt = 
5MHz are used, a, Homodyne signal generated by the out- 
put beam as the input power is slowly swept up and down 
over the Hopf bifurcation point pi. The gray bars labeled 
with 'b' and 'c' indicate the regions of which the zoomed-in 
views are displayed in b and c. d, Theoretical prediction of 
the input-output characteristics of the system. The solid and 
the dotted black lines represent the stable and the unstable 
steady state solutions, respectively. Points pi and p2 indicate 
Hopf bifurcations discussed in the text. The upper and lower 
blue lines trace the maximum and the minimum output power 
of the limit cycles observed in simulations, e, On the bottom, 
the predicted output power of the system (from d) near the 
subcritical bifurcation point pi is compared with the mea- 
surement. The measured data is plotted in red (up-sweep) 
and green (down-sweep), using the following rules: When no 
oscillation is detected, the mean output power is recorded by 
a cross-mark. When the oscillation is detected, a sinusoid is 
fitted to the measured signal and its max/min optical power is 
denoted by a pair of connected circles. On the top, we plot the 
frequencies of the both simulated and measured limit cycles. 
The simulation result is plotted in brown, while the measure- 
ments are denoted by red (up-sweep) and green (down-sweep) 
square marks. 



limit cycle. With the system oscillating, the input power 
can be lowered back down and the limit cycle remains 
until another bifurcation occurs at a lower input power 
(^1150nW) than pi. As shown in Fig. [2^, such hystere- 
sis can be observed clearly in our experiment, and we 
find that the measured amplitude and frequency of the 
limit cycle closely match predictions. Moreover, com- 
bined measurements of the amplitude and phase quadra- 
tures of the output field confirm that the oscillation of 



the intracavity field follows a trajectory on the optical 
phase plane that is well predicted by theory (See Supple- 
mentary Fig. SI). 

Near the predicted super-critical Hopf bifurcation 
point p2 (~2100nW), on the other hand, the match be- 
tween our measurements and theory was not as good. 
Simulation predicts that the amplitude of the limit cy- 
cle should gradually converge to zero as the system ap- 
proaches p2 from below, at which point a stable equilib- 
rium point defines the steady-state. However, the experi- 
mental system continued to show large-amplitude oscilla- 
tory behavior at random intervals in time even when the 
input power was far above the predicted location of p2. 
Simulations with added technical noise suggest that such 
behavior is quite reasonable, however, as the equilibrium 
points near p2 are highly under-damped. 

The limit cycle behavior of our system converts the 
constant input power into an rf-modulated output power. 
We are thus able to demonstrate an elementary physical 
mechanism that potentially could be adopted to make 
self-contained inline optical oscillators in a nanopho- 
tonic setting, for which the use of conventional optical- 
oscillator-type devices such as Q-switched lasers might 
not be practical. Our system furthermore is capable of 
transferring baseband power modulation to rf modula- 
tion of the output via the correlation between the in- 
put power and the frequency of the limit cycle (Fig. [2^, 
upper plot). This suggests the possibility of an all- 
optical baseband-to-rf up-converter. And finally, hystere- 
sis around the sub-critical Hopf bifurcation can be used 
to make a latch of rather unique characteristics, which 
switches between dc and modulated signal formats. Op- 
tical dc-rf latching devices may find use in optical com- 
munication systems where signal up-conversion is a com- 
mon practice to avoid the high level of the background 
noise at low frequencies. Fig. [3] shows a demonstration of 
such latching performance. In this context we note that 
solid-state implementations of this type of cavity nonlin- 
ear optical device, for example using ensembles of quan- 
tum dots [33] or vacancy centers [34] , could achieve much 
higher values of (g , 7_l , k) and therefore much higher 
limit cycle frequencies, potentially in the GHz range. 

We find that the energy stored in the optical dc-rf 
latch in Figj3] is comparable to that of the elementary 
latch in Fig. ^ as in both cases the cavity contains about 
500-1000 photons (100-200 attojoule) on average in the 
high energy state and much fewer when in the low en- 
ergy state. However, the energy required for switching 
between the two states was more than an order of mag- 
nitude higher for the dc-rf latch (~30 fJ). This is mostly 
due to the mismatch between the input laser frequency 
and the resonances of the system. In Supplementary Sec- 
tion Fig. S2 we show a demonstration of sub-fcmtojoule 
dc-to-rf switching operation of the same Hopf bifurca- 
tion latch using a secondary near-resonant beam. We 
also note that the use of pulse shaping on the control 
pulse |35j could potentially reduce the switching energy 
of both the elementary and the dc-rf latch further. 
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FIG. 3. Demonstration of a dc-rf latch. The system 
shown in Fig.[2]is used to demonstrate latching between a sta- 
ble equilibrium point and a limit cycle. The spectrogram in 
the middle is generated using 0.2/is segments of the homodyne 
data, which contain roughly 5 oscillation cycles. OdB power 
level was arbitrary chosen. Note that more bandwidth is re- 
tained in the homodyne signal displayed here than in Fig.[]J, 
in order to preserve the limit cycle oscillations. In the bottom 
frame we show a zoomed-in view of the homodyne signal from 
3/is to 4/xs. 



In conclusion, we have demonstrated that a cavity- 
based optical device that incorporates a resonant non- 
linear response can exhibit complex dynamic behaviors 
such as hysteretic limit cycle formation at femtojoule en- 
ergy scales. We have suggested novel signal processing 
functions based on such phenomena, and have demon- 
strated an all-optical configuration that simultaneously 



latches a logic level and performs rf modulation of the 
output power. 
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I. METHODS 

The atoms are first collected and cooled by a magneto- 
optical trap (MOT) formed on a gold-coated mirror sur- 
face [36]. This surface MOT is located outside the cav- 
ity, 2cm away from the cavity axis, where the cooling 
beams are not obstructed by the cavity mirrors. The 
atoms are then further cooled to sub-Doppler temper- 
ature (< 10/iK) by polarization-gradient cooling, opti- 
cally pumped to the mj?=4 state, and then loaded in the 
magnetic trap. The magnetic trap, formed by currents 
through copper wires (100 — 200/im diameter) buried un- 
derneath the mirror surface, can be moved by shifting 
the currents [23] . and is used as a conveyer to trans- 
port atoms (~ 10 6 ) into the cavity mode. At the end of 
the transport process (200 — 300ms), the atoms are re- 
leased from the magnetic trap and a large uniform mag- 
netic field is applied in the direction of the cavity axis, 
which ensures that the quantization axes for the atoms 
and the field are aligned in the presence of possible stray 
fields near the surface. The strong bias field also lifts 
the degeneracy among the Zeeman states (Zeeman shift: 
2ir ■ 17MHz). Once the atoms are in position, the input 
laser is injected from one side of the Fabry-Parot cavity 
and drives the TEM o mode. The beam that leaks out 
from the other side of the cavity serves as the output of 
the system. Due to the loss at the cavity mirrors and the 
existence of two possible directions for the leakage, we es- 
timate rj, the probability of an intracavity photon making 
it to the output channel, to be 0.41. The output beam 
goes through an optical isolator and is then measured by 
a homodyne/heterodyne detection scheme fTU]. 
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